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The origin of the G-type antiferromagnetism (AFM(G)) and puzzling properties of RTiOs with 
_R=La are studied. We clarify that the crystal field from La caused by the GdFe03-type distortion 
lifts the t2 g degeneracy at Ti 3d orbitals. The lift stabilizes the AFM(G) with spin-exchange constant 
in agreement with neutron scattering results. The orbital-spin structures for R=Pr, Nd and Sm 
are also consistent with experiments. We propose that the GdFeOs-type distortion has a universal 
mechanism of controlling orbital-spin structure competing with the Jahn- Teller (JT) mechanism. 



Rich magnetic and orbital phases, their phase tran- 
sitions and interplays have attracted much interest in 
strongly correlated electron systems, particularly in tran- 
sition metal (TM) oxides. Coupling to lattice further en- 
riches the interplay through lattice distortions, dynam- 
ical phonons and cooperative effects such as JT distor- 
tions 0. Cuprate superconductors and manganese per- 
ovskite compounds with colossal magnetoresistance be- 
long to the TM oxides with 3d e g bands at the Fermi 
level. Perovskite-type Ti oxides i?TiC>3 with R being 
a rare-earth ion belong to the same type but with 3d 
tig bands at the Fermi level and show a very different 
behavior. Particularly, while the above mother materi- 
als are commonly typical Mott insulators, i?TiC>3 shows 
markedly rich and complicated behaviors possibly due to 
the above interplay. The crystal structure of i?Ti03 is a 
pseudocubic perovskite with an orthorhombic distortion 
(GdFe03-type distortion) and exhibits a magnetic phase 
transition as a function of the magnitude of this distor- 
tion HUSQ. The end compounds R =Gd and Y with 
large GdFe03-type distortions showing the ferromagnetic 
(FM) ground state are relatively well understood |(J • 
In contrast, the magnetic structure of LaTiC-3 with rela- 
tively small GdFe03-type distortion has AFM(G) struc- 
ture with Neel temperature (T N ) of - 130 K Q. The 
origin of this G-type antiferromagnetism has long been 
puzzling and controversial. The purpose of the present 
theoretical study is to solve this puzzle and offer an over- 
all understanding of the interplay in these compounds. 

In perovskite TM compounds, the JT distortions often 
play crucial roles in determining low-energy electronic 
states. However, with early TM ions, the JT coupling is 
considerably weaker than that in the late-3c? compounds. 
Actually, the distortion of the Ti06 octahedra in LaTi03 
has not been detected in contrast to the manganites and 
cuprates Therefore, in LaTi03, the crystal field of 

O ions surrounding the Ti 3+ ion has a cubic symmetry so 
that the degeneracy of the ti g level survives. Under this 
circumstance, the AFM(G) ground state is surprising and 
controversial since in the orbitally degenerate system, we 
expect that a FM state with antiferro-orbital ordering is 
stabilized both by transfers and the Coulomb exchange 
interaction yj. Indeed, a recent weak coupling study 
shows that a FM state, out of which two states (yz + 
izx) I V2 | and xy f are alternating is favored both by the 
relativistic spin-orbit (LS) interaction and by the spin- 



orbital superexchange interaction |l0j . 

Moreover, a recent neutron-scattering study shows the 
spin-wave spectrum well described by an isotropic spin- 
1/2 Hciscnbcrg model with a nearest-neighbor superex- 
change constant J ~ 15.5 meV .1.1]. If the orbital mo- 
ment is nonzero, the anisotropy appears through the LS 
interaction. Therefore, the isotropy indicates that the 
orbital moments are quenched. 

Recently, a possible orbital liquid state was proposed 
by assuming the G-type antiferromagnetism on the ba- 
sis of resultant small exchange interaction in the orbital 
sector y\^. However, the origin of the AFM(G) state in 
LaTi03 is still puzzling when the orbital is disordered. 
On the other hand, the possibility of the D 3( i distortion 
of the TiOe octahedra was also proposed for the origin 
of the AFM(G) 13]. In this distortion, the Ti0 6 octahe- 
dron is contracted along the threefold direction, and the 
threefold-degenerate i2g-levels split into a nondegenerate 
lower aig-level and twofold-degenerate higher e g -levels. 
The lift of the ti g levels with the occupation of the lowest 
ai g -orbital well explains the emergence of the AFM(G) 
and the isotropic spin- wave spectrum. However, this D^d 
distortion has so far not been detected. 

In addition to the JT distortion, the GdFe03-type lat- 
tice distortion is another generic phenomenon in per- 
ovskites. It has been assumed mainly to control the band- 
width through the M-O-M angle with M being TM ions, 
while its direct effects on interplay of spins and orbitals 
have not been considered seriously. 

In this letter, by taking a particular example of 
LaTi03, we clarify that the generic GdFe03-type dis- 
tortion actually generates a new mechanism for control 
of orbital-spin low-energy structure through lift of the or- 
bital degeneracy by crystal fields of R ions. This mech- 
anism competes with that of the JT distortions and LS 
interactions. By utilizing the experimentally obtained 
coordination parameters, we construct the Hamiltonian 
for the crystal field of La cations in LaTiOs. The anal- 
yses of the obtained Hamiltonian show that the shifts 
of La cations due to the GdFe03-type distortion turn 
out to generate the crystal field which is similar to the 
D3d-crystal-field. As a result, the threefold degeneracy of 
the cubic-t2g-levels splits into three nondegenerate levels. 
The calculations of the energies and the spin-exchange 
constant well explain the emergence and the properties 
of the AFM(G) state in LaTi03. In addition, we also 
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FIG. 1: (a) The distortions of the R cations along the b- 
axis due to the GdFeC>3-type distortion are presented. In 
_RO-plane 1, they distort in the positive direction while in 
_RO-plane 2, they distort in the negative direction, (b) In the 
case of site 1, the distances between the Ti ion and R ions 
located in the ±(1, 1, — Indirections are decreased. As a re- 
sult, the La ions generate the crystal field which is similar 
to the £>3d-crystal field with [1, 1, — l]-trigonal axis (see text). 
The distortions of the O ions along the c-axis which are in- 
duced by the tilting of the octahedron are also presented, (c) 
The energy-level structure in the crystal field due to the R 
cations. The threefold degeneracy of the cubic-t2 9 levels split 
into three nondegenerate levels. 

study the effects of the crystal field of rare-earth ions in 
i?Ti03 with R being Pr, Nd and Sm. Again experimen- 
tally observed orbital structure as well as the reduction 
of spin exchange are reproduced. We point out the im- 
portance of GdFeC>3-type distortion as a generic control 
mechanism of orbital-spin structure in perovskite com- 
pounds. 

In the GdFe03-type distortion, the R ions distort 
mainly along the (1, l,0)-axis or the 6-axis, and slightly 
along the (1, —1, 0)-axis or the a-axis. In Fig. ^ (a), we 
show the distortions along the 6-axis. There are two 
kinds of i?0-planes (plane 1 and plane 2) stacking alter- 
natingly along the c-axis. In plane 1, the R ions shift 
in the positive direction while they shift in the nega- 
tive direction in plane 2. Consequently, the crystal field 
from the R cations are distorted from a cubic symme- 
try. For example, in the case of sites 1 and 2, the dis- 
tances between the Ti ion and the R ions located in 



the ±(1,1, — Indirections decrease while those along the 
±(1, 1, Indirections increase (see also Fig.^(b)). On the 
other hand, the distances along ±(1,1, Indirections de- 
crease while those along ±(1, 1, — l)-directions increase 
for sites 3 and 4. The changes of the other Ti-i? dis- 
tances are rather small. Since the nominal valence of R 
ion is 3±, we expect that Ti 3d orbitals directed along 
the shorter Ti-R bonds are lowered in energy because of 
the attractive Coulomb interaction. Without any distor- 
tions of the Ti06 octahedra, the crystal field of the ligand 
oxygens has a cubic symmetry. When we introduce the 
crystal field of R cations (Hm) as a perturbation, the 
threefold degeneracy of the cubic-i2gTevels split. Here, 
for simplicity, we take Hm by assuming that there ex- 
ist the point charges with +3 valence on each R cation. 
The Coulomb interaction between an electron on a Ti 3d 
orbital and a La 3+ ion is given by using the dielectric 
constant £Tiii as, 



where expresses the coordinates of the i-th R ion, 
and Zr(=+3) is the nominal valence of La +3 ion. We 
calculate the following matrix elements, 

(m\v(v)\m!) = J drtp* 3dm v(r)ip 3d , m > (2) 

with 

¥>3d,m = R3d(r)Y 2m {0,<l)). (3) 

Here, the integer m(=— 2, ...,2) is the magnetic quantum 
number. The coordinates of the La ions are calculated 
by using the positional parameters and the cell constants 
obtained by the x-ray diffraction study @. After trans- 
forming the basis from m to the t^g representations, we 
obtain the Hamiltonian Hm with xy-, yz- and zx-basis. 
These representations are defined in the coordinates with 
x-, y- and z-axes attached to each TiOe octahedron. 

By diagonalizing the obtained Hm, we can evalu- 
ate the 3d-\eve\ energies and their representations. The 
threefold degeneracy of the t^g orbitals is split into three 
isolated levels as shown in Fig. ^ ( c ) ■ The energy lev- 
els Ai and A 2 are 0.7685/£ T iLa eV and 1.5849/e T1La 
eV, respectively. The representations of the lowest lev- 
els at each site can be specified by the linear combi- 
nations of xy, yz and zx orbitals as axy — byz — czx, 
axy — cyz — bzx, axy + byz + czx and axy + cyz + bzx, re- 
spectively, where a 2 +fe 2 +c 2 =l. Here, we note that there 
exists a mirror plane vertical to the c-axis, and the orbital 
structure has the same symmetry as the orthorhombic 
GdFe03-type distortion. The coefficients take the values, 
a=0.6046, 6=0.3900 and c=0.6945. This orbital structure 
is similar to that in the Dsd-crystal-fields with [1,1,-1]- 
, [1,1,-1]-, [1,1,1]- and [l,l,l]-trigonal axes at sites 1, 
2, 3 and 4, respectively. The representations of the low- 
est levels in these Z^-crystal- fields are (xy — yz — zx) , 

A=(xy-yz-zx), -j=(xy + yz + zx) and -j=(xy + yz + zx), 
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respectively. In the previous study 0, this configura- 
tion was referred to as config. 4, and it was shown that 
the D3d-crystal-field with this trigonal-axes configuration 
strongly stabilizes the AFM(G) spin structure. 

We next discuss the stability of the magnetic state 
in the crystal field of Hr\. For this purpose, we em- 
ploy the effective spin and pseudospin Hamiltonian Q: 
H c s. — -ff cry . + H ST . The first term -ff cry . denotes the 
crystal field on the Ti ti g orbitals. The second term 
H ST is the spin and pseudospin term, in which the three- 
fold tig orbital degrees of freedom are represented by 
the pseudo-spin-1 operator r. This effective Hamiltonian 
is derived from the multiband d-p model in the insulat- 
ing limit by following an approach similar to the Kugel- 
Khomskii formulation [l4lll5l |. In this d-p model, the full 
degeneracies of Ti 3d and O 2p orbitals as well as the on- 
site Coulomb and exchange interactions are taken into 
account Q. The on-site Coulomb interactions are ex- 
pressed using Kanamori parameters u, y! , j and j' which 
satisfy the following relations 0, 0: u = U + ^-j, 
vf = u — 2j and j = j' . Here, U gives the magnitude 
of the multiplet-averaged d-d Coulomb interaction. The 
charge-transfer energy A, which describes the energy dif- 
ference between occupied O 2p and unoccupied Ti 3d 
levels, is defined using U and the energies of the bare 
Ti 3d and O 2p orbitals e d and e p as A = e d + U — e p , 
since the characteristic unoccupied 3d level energy on the 
singly occupied Ti site is e d + U. The values of the pa- 
rameters are estimated by the cluster-model analyses of 
valence-band and transition-metal 2p core-level photoe- 
mission spectra and the analyses of the first-principles 
band calculations [l8L [Tsj . We take the values of these 
parameters as A = 5.5 eV, U = 4.0 eV, V p da = -2.4 eV, 
VpdK = 1.3 eV, V PP<J = 0.52 eV, V pp7T = -0.11 eV and 
j = 0.46 eV. The effects of the CdFe0 3 -type distortion 
arc considered through the transfer integrals which are 
defined using Slater-Koster parameters [2(|. Substitut- 
ing the -ffcry. term with -Hri, we calculate the energies for 
several magnetic structures by applying the Hartree-Fock 
approximation. In Fig. [2 we plot the calculated energies 
as functions of the magnitude of Ai (=0.7685/etlr)- We 
tune the magnitude of Ai by varying Stir- In the region 
of Ai >0.03 eV, the AFM(G) state is strongly stabilized 
relative to the other structures. In this region, Ai is 
much larger than k^T^ so that the orbital occupation 
is restricted to the lowest level irrespective of the spin 
structure. We can estimate the magnitudes of the spin- 
exchange interactions for each Ti-Ti bond in the subspace 
of singly-occupied lowest levels. The spin-exchange con- 
stant J can be represented as J = (E^ — E^i)/2S 2 with 
an d ^Tl being the energy gains of the Ti-Ti bond 
for ||- and fl-pairs, respectively. For LaTi03, the val- 
ues of J along the x-, y- and z-axes take as =18.45 
meV, J y =18.45 meV and J z = 19.71 meV, respectively. 
They are in agreement with the value obtained by the 
neutron scattering experiment of ~15.5 meV, and are 
consistent with the spin-wave spectrum well described 
by the isotropic Heisenberg model. 

We further note that the crystal field of the La cations 
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FIG. 2: Energies for several magnetic structures under the 
crystal field of Hri are plotted as functions of Ai. 



has two origins. One is the Coulomb potential from 
charged La ions as we have studied above, and the other 
is the hybridization between the Ti 3d orbitals and un- 
occupied orbitals on the La cations. We next estimate 
the effect of the crystal field due to the latter origin for 
LaTiOs. We construct the Hamiltonian for hybridization 
between Ti 3d and La 5d orbitals (Hri) by using the 
second-order perturbational expansion in terms of the 
transfers between Ti 3d and La 5d orbitals (t dd ). The 
expression of the matrix element of EIri is 

j.dd j.dd 

{m\H m \m') = -Y J m "' ; ' 7 (4) 

Here, the indices m and m! run over the cubic-^g- 
representations, xy, yz and zx. The symbols i and 7 
are indices for the eight nearest-neighbor La ions and the 
fivefold La 5d orbitals, respectively. The symbol A$d de- 
notes the characteristic energy-difference between Ti tig 
and La 5d orbitals. The transfers t dd are given in terms 
of Slater-Koster parameters Vdda, Vdd-K and VddS- It is 
assumed that these parameters are proportional to d~ 5 
with d being the Ti-La bond length. On the basis of the 
analyses of LDA band structure, we fix these parameters 
as V dda = -1.04 eV, Vmk = 0.56 eV and V dd s = eV 
for Ti-La bond length of 3.5 A, and A 5d = 3.6 eV [Hj. 
By diagonalizing thus obtained Hm, we have evaluated 
the energy- levels and their representations (see Table QJ. 
In the obtained energy-level structure, the threefold tig- 
levels split into three nondegenerate levels, which is simi- 
lar to that of Tdfji . Since the representation of the lowest 
level of Hjii is also similar to that of Hri , the level split- 
ting of Hri+Hri is well expressed by the sum of these 
two contributions. Considering the fact that AFM(G) 
state is stabilized in the region of Ai >0.03 eV, and 
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TABLE I: The energy-level structures and the representations 
of the lowest levels for Hri and Hr2- 







Hri 


Hr2 




Hri 


Hr2 


Ai 


(eV) : 


0.7685/ETiLa 


0.0849 


a : 


0.6046 


0.6233 


A 2 


(eV) : 


1.5849/era* 


0.1443 


b : 


0.3900 


0.4385 










c : 


0.6945 


0.6474 



Ai for Hb2 is ^0.085 eV, the crystal field due to the 
hybridization between Ti 3c? and La 5d orbitals (Hrs) 
alone turns out to stabilize the AFM(G) spin structure 
strongly. In addition, the value of Ai is sufficiently large 
as compared with the coupling constant of the LS in- 
teraction in Ti 3+ (Cd=0.018 eV) 22\. Consequently, the 
crystal field of La cations dominates over the LS inter- 
action, resulting in the quenched orbital moment. The 
reduction of the ordered magnetic moment is attributed 
to the itinerant fluctuation instead of the LS interaction 
near the metal-insulator phase boundary as discussed in 
the literature [lflj . 

We have also studied the effects of the crystal field of R 
ions in i?Ti03 with R being Pr, Nd and Sm. The analy- 
ses of the crystal field Hamiltonian again reproduces the 
AFM(G) ground state in each compound. Moreover, the 
gradual decrease of the spin-exchange constant J as R 
goes from La, Pr, Nd to Sm is reproduced (for example J 
shows 20% reduction at Nd from La), which is consistent 
with the gradual decrease of T/v of these compounds. The 
orbital states show that the occupation of the yz orbital 
decreases at sites 1 and 3 as the GdFeOvtype distortion 
increases while at sites 2 and 4, the zx-occupation de- 
creases. As a result, the orbital structure in which sites 
1, 2, 3 and 4 are occupied by -^(xy — zx), -^{xy — yz), 

-^(xy + zx) and -^(xy + yz), respectively is progres- 
sively stabilized. For example, for Sm, the coefficients 
take a = 0.63, b = 0.73 and c = 0.24 even without the JT 
distortion. It has the similar symmetry as that realized 



in YTiC>3 with a large JT distortion, which is consistent 
with recent RXS study j2^|. It is interesting to note that, 
in our results, the character of the orbital symmetry in 
YTiC>3 is retained but quantitatively and continuously 
decreases toward the type of LaTiC>3 with the decreasing 
GdFe03-type distortion. When the crystal field from O 
ion due to the JT distortion is superimposed, our calcu- 
lation shows that the dominance is taken over by the JT 
mechanism at R =Sm. Through severe competition at 
R =Sm, the Jahn- Teller mechanism comes to control the 
orbital-spin structure at R —Gd and Y. 

In summary, in order to elucidate the origin of the 
AFM(G) state in RTi0 3 with R=Ti, we have examined 
the crystal field of La cations by considering the exper- 
imentally observed distortion of the La ions. The crys- 
tal field Hamiltonian is constructed by using the experi- 
mental position parameters. Based on this Hamiltonian, 
we have shown that the distortion of the La ions caused 
by the GdFe03-type distortion generate the crystal field 
with lower symmetry, which is similar to the D^d crystal 
field. Then the threefold degenerate cubic-i2 ff -levels split 
into three isolated levels. The energies and the spin- 
exchange constant calculated by the effective Hamilto- 
nian show that the lowest-level occupations in this crys- 
tal field stabilize the AFM(G) state, and account for the 
spin-wave spectrum for LaTi03 with (1) J ~ 15.5 meV, 
(2) isotropic spin-coupling, and (3) quenched orbital mo- 
ment. The orbital-spin structre for R = Pr, Nd and Sm is 
also accounted for by the same mechanism. GdFe03-type 
distortion has a universal relevance in determining the 
orbital-spin structure of perovskite compounds in com- 
petition with the JT mechanism. 

We thank N. Hamada for providing us with data 
on LDA calculations of LaTiOs. This work is sup- 
ported by "Research for the Future Program" (JSPS- 
RFTF97P01103) from the Japan Society for the Promo- 
tion of Science. 



[1] For a review see M. Imada, A. Fujimori and Y. Tokura: 

Rev. Mod. Phys. 70 1039 (1998). 
[2] Y. Okimoto, T. Katsufuji, Y. Okada, T. Arima and Y. 

Tokura: Phys. Rev. B 51 9581 (1995). 
[3] T. Katsufuji, Y. Taguchi and Y. Tokura: Phys. Rev. B 

56 10145 (1997). 
[4] J. P. Goral, J. E. Greedan and D. A. MacLean: J. Solid 

State Chem. 43 244 (1982). 
[5] J. E. Greedan: J. Less-Common Met. Ill 335 (1985). 
[6] M. Mochizuki and M. Imada: J. Phys. Soc. Jpn. 69 1982 

(2000) . 

[7] M. Mochizuki and M. Imada: J. Phys. Soc. Jpn. 70 1777 

(2001) . 

[8] D. A. MacLean, H.-N. Ng and J. E. Greedan: J. Solid 

State Chem. 30 35 (1979). 
[9] M. Eitel and J. E. Greedan: J. Less-Common Met. 116 

95 (1986). 

[10] M. Mochizuki: J. Phys. Soc. Jpn. 71 2039 (2002). 



[11] B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M. v. Zim- 
mermann, J. P. Hill, D. Gibbs, Y. Taguchi and Y. Tokura: 
Phys. Rev. Lett. 85 3946 (2000). 

[12] G. Khaliullin and S. Maekawa: Phys. Rev. Lett. 85 3950 

(2000) . 

[13] M. Mochizuki and M. Imada: J. Phys. Soc. Jpn. 70 2872 

(2001) . 

[14] K. I. Kugel and D. I. Khomskii: Pisma. Zh. Eksp. Teor. 

Fiz. 15 629 (1972). 
[15] K. I. Kugel and D. I. Khomskii: Zh. Eksp. Teor. Fiz. 64 

1429 (1973). 
[16] B. H. Brandow: Adv. Phys. 26 651 (1977). 
[17] J. Kanamori: Prog. Theor. Phys. 30 275 (1963). 
[18] T. Saitoh, A. E. Bocquet, T. Mizokawa and A. Fujimori: 

Phys. Rev. B 52 7934 (1995). 
[19] P. Mahadevan, N. Shanthi and D. D. Sarma: Phys. Rev. 

B 54 11199 (1996). 
[20] W. A. Harrison: Electronic Structure and the Properties 



■5 



of Solids (Dover, New York, 1989). 
[21] N. Hamada (private communication) 
[22] S. Sugano, Y. Tanabe and H. Kamimura: Multiplets of 

Transition-Metal Ions in Crystals (Academic, New York 

1970) 



[23] M. Kubota, H. Nakao, Y. Murakami, Y. Taguchi, Y. 
Tokura, M. v. Zimmermann, J. P. Hill and D. Gibbs: 
unpublished. 



